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We prove that the oft-used stationary-phase method gives a very accurate expression for the Fourier 
transform of the gravitational-wave signal produced by an inspiraling compact binary. We give 
three arguments. First, we analytically calculate the next-order correction to the stationary-phase 
approximation, and show that it is small. This calculation is essentially an application of the 
steepest-descent method to evaluate integrals. Second, we numerically compare the stationary- 
phase expression to the results obtained by Fast Fourier Transform. We show that the differences 
can be fully attributed to the windowing of the time series, and that they have nothing to do with 
an intrinsic failure of the stationary-phase method. And third, we show that these differences are 
negligible for the practical application of matched filtering. 
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I. INTRODUCTION AND SUMMARY 

The ongoing construction of kilometer-size interfero- 
metric gravitational-wave detectors, such as the Ameri- 
can LIGO (Laser Interferometer Gravitational-wave Ob- 
servatory), the French-Italian VIRGO, and the British- 
German GEO600, has motivated a lot of recent work on 
strategies to analyze the data to search for and measure 
the gravitational- wave signals |H-|l6f| . Much of this work 
has been devoted to inspiraling compact binaries, com- 
posed of neutron stars and/or black holes, which are one 
of the most promising sources of gravitational waves for 
these detectors. The key idea behind this work is that 
the gravitational-wave signals will be detected and mea- 
sured by matched filtering [jlTj, a well-known technique 
by which the detector's data stream is cross-correlated 
with a set of model waveforms (called templates), and 
the signal-to-noise ratio maximized as a function of the 
template parameters. 

While the detector output is represented by a discrete 
time series, the operations associated with matched filter- 
ing are usually carried out in (discrete) frequency space, 
which requires taking the (discrete) Fourier transform of 
the time series; the method of choice here is the standard 
Fast Fourier Transform (FFT). At the same time, it is 
necessary to compute the Fourier transform of h(t), the 
model waveform. The waveform h(t) can easily be dis- 
cretized to mimic the discrete sampling of the detector 
output, and the discrete time series can easily be Fourier 
transformed by FFT to yield the discrete analogue of the 
(continuous) Fourier transform h(f). 

However, in theoretical investigations [p]- ^3|Jl^Jlr| it is 
often much more convenient to deal with h(f) as an an- 



alytic expression rather than as a table of values; efforts 
to obtain such an analytical expression, even if it is only 
approximate, are therefore well justified. It can also be 
argued that to express the waveform in the frequency do- 
main is in a sense more natural than to express it in the 
time domain. One of the reasons is that the orbital en- 
ergy, whose expression is required in the derivation of the 
waveform, is primarily a function of F, the gravitational 
wave's instantaneous frequency (which will defined pre- 
cisely below) ; it is therefore natural to express h also as a 
function of F, which unlike t is a coordinate-independent 
quantity. Another reason resides with the fact that the 
post-Newtonian expression for the relation F{t) suggests 
that the instantaneous frequency is not always a mono- 
tonically increasing function of time jl8| , a nonphysical 
behavior that signals the eventual breakdown of the post- 
Newtonian expansion; in contrast, the post-Newtonian 
expression for the inverted relation t(F) is monotonic, a 
property that suggests that F is indeed a more natural 
time variable. 

In the past, and in the specific context of inspiral- 
ing binaries, analytic expressions for h(f) have been ob- 
tained within the stationary-phase approximation [fl9|| , 
essentially the leading-order term in an expansion in 
powers of the small quantity (radiation-reaction time 
scale) /(orbital period). While it has generally been felt 
that this approximation is adequate, this belief has not 
yet been backed up (in the published literature) with a 
detailed quantitative analysis. Claims in the literature 
|l4| to the effect that this approximation is not adequate 
have prompted us to examine this question. 

Our objective in this paper is to prove, beyond any rea- 
sonable doubt, that the stationary-phase approximation 
to the Fourier transform of an inspiraling-binary signal 
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is in fact very accurate. First, we calculate the next- 
order correction to the stationary-phase approximation, 
and show that it is indeed small. This calculation is es- 
sentially an application of the steepest-descent method 
to evaluate integrals (l9). We then numerically compare 
the stationary-phase expression for h(f ) to the results 
obtained by FFT. We show that the differences can be 
fully attributed to the windowing of the time series, and 
that they are irrelevant for matched filtering. 

We now present and explain our main results. As we 
shall justify below, it is sufficient for our purposes to con- 
sider the leading-order expression for the gravitational- 
wave signal, obtained by assuming that the binary's or- 
bital motion (taken to be circular) is governed by New- 
tonian gravity, with an inspiral — caused by the loss 
of energy and angular momentum to gravitational radi- 
ation — governed by the Einstein quadrupolc formula. 
This so-called Newtonian signal is given by jjj 



M 



h(t) = Q(angles) — (ttMF) 2 / 3 e 
r 



(1.1) 



where Q is a function of all the relevant angles (posi- 
tion of the source in the sky, inclination of the orbital 
plane, orientation of the gravitational- wave detector), 
M = (mim2) 3 ' 5 /( m i + rn-if-l^ (with m\ and 1712 de- 
noting the individual masses) is the chirp mass, r is the 
distance to the source, F(t) the instantaneous frequency 
(twice the orbital frequency), and $(i) = J 2nF(t') dt' is 
the phase. The function hit) represents the gravitational 
wave measured at the detector site. It is a linear superpo- 
sition of the wave's two fundamental polarizations, and 
we choose to express it in this complex form for conve- 
nience. The relation between F and t is given by |]l| 
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where t c , the "time at coalescence", is such that formally, 
F(t c ) = 00. The relation between $ and t is 
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where $ c , the "phase at coalescence", is equal to &(t c )- 
The Newtonian signal is the leading-order term in the 
expansion of the gravitational waves in powers of V <C 1 , 
where V — (nMF) 1 / 3 is (up to a numerical factor) the 
orbital velocity. Post-Newtonian corrections to this re- 
sult come with a relative factor of order V in the ampli- 
tude, and a relative factor of order V 2 in the phase p0[ . 
Throughout this paper we use geometrized units, setting 
G = c= 1. 

In the stationary-phase approximation, the Fourier 



transform of the function h(t) appearing in Eq. (1.1) is 
given by || 



h spa .(f) = 



/ 30tt QM 2 -t/2 
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where v = (irMf) 1 ^ 3 <C 1, and 



ip(v) = 2nft c -$ c - - 
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(1.5) 



Two types of corrections to this result are calculated in 
Sec. II. 

The first type of correction constitutes an intrinsic im- 
provement on the stationary-phase approximation. We 
show that the steepest-descent evaluation of the Fourier 
transform returns the same amplitude as before, but that 
the phase is altered by a term Sip: ip — > %j) + Sip, where 
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6ip(v) = l- v 5 + 0(v 10 ). 
45 



(1.6) 



Notice that this is a correction of order v relative to 
ip(v). This is much smaller than post-Newtonian correc- 
tions to the phase, which appear at relative order v 2 [pj. 
Incorporating these post-Newtonia n co rrections into our 
calculation would only change Eq. (1.6) by adding a term 
of order v 7 to the right-hand side. This justifies the fact 
that it was not necessary, for the purposes of t his i nvcs- 
tigation, to use more accurate versions of Eqs. (1.2) and 

The second type of correction addresses an implicit 
assumption of the stationary-phase method, that the 
function h(t) has support in the complete time interval 
—00 < t < t c . Physically, this assumption means that 
the binary system must have formed in the infinite past 
(a reasonable assumption given the long lives of com- 
pact binaries), and that the inspiral must continue until 
F = F(t c ) = 00 (an unrealistic assumption). If we choose 
instead to restrict the time interval to i m i n < t < < max , 
such that F min = F(t mhl ) > and F max = F(t max ) < 00, 
then the Fourier transform will be affected, and it will 
differ from h spa (f). The value of F m - ln is typically chosen 
to reflect the lower bound of the instrument's frequency 
band. The value of F max could be chosen to reflect the 
upper bound of the instrument's frequency band, or for 
some binaries it can be chosen to correspond to the ap- 
proximate frequency of the last stable orbit, at which the 
inspiral signal changes over to a poorly-known merger 
signal. 

We shall refer to this truncation of the time interval 
as "windowing" , and for concreteness, we will assume 
that the signal is started abruptly at t — £ m i n and ended 
abruptl y at t — t niax . Thus, h(t) is assumed to be given 
by Eq. (1.1) in the interval i m ; n < t < < max , and is as- 
sumed to be zero outside this interval. In Sec. II we 
show that windowing affects both the amplitude and the 
phase of the Fourier transform. The amplitude acquires 
an extra factor 1 + 8 A„ , where 
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where x m i n = (^MF^^ 1 ^ 3 and 
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with similar equations holding for x max and </> max . On 
the other hand, the phase acquires an extra term 5?p w 
given by 



Sipw{v) 
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The fact that SA W and both diverge at the bound- 
aries V — — Xniin 

and v = s max signals the breakdown of our 
approximations there. Away from the boundaries, <5A W 
and 5tp w are bounded, and they oscillate as a function of 
frequency. Thus, these corrections represent amplitude 
and phase modulations induced by the abrupt cutoffs of 
the function h(t) at the boundary points. This is to be 
contrasted with Sip, which represents a steady phase drift. 

Our complete expression for the steepest-descent ap- 
proximation to the Fourier transform is therefore 



W/) =U/)( 1 + M,)e 



i(5ip-\-8ip v . 



(1.10) 



In Sec. Ill we show that the only noticeable corrections 
to h spa (f) are the amplitude and phase modulations that 
come as a consequence of windowing; in particular, the 
intrinsic correction Sip is too small to be noticeable in 
the relevant frequency interval. We do this by comparing 
hap&(f) to hfft(f), the discrete Fourier transform of the 
windowed time series h(t); this comparison reveals that 
any discrepancy between the two versions of the Fourier 
transform can be fully accounted for by the modulations 
6A W and 5ip w . This allows us to conclude that window- 
ing, and windowing only, must be held responsible for any 
discrepancy between /i spa (/) and /ifft(/)- While limited 
to the Newtonian signal of Eq. ([0]), there is no reason 
to believe that this conclusion would be invalidated by a 
full post-Newtonian analysis pl[ . 

Finally, in Sec. IV we calculate the matched-filtering 
overlap between h spa (f) and /ifft(/), and show that the 
modulations do not significantly affect the overlap. This 
result, together with our previous findings, lead us to con- 
clude that for the purposes of matched filtering, h spa (f) 
and /ifft(/) are essentially equivalent representations of 
the gravitational-wave signal. 

The main conclusion of this work is that the stationary- 
phase method returns a sufficiently accurate expression 
for h(f). What's more, from the fact that Sip/ip — 
0(v 10 ), we can be sure that the method will stay ac- 
curate for as long as the post-Newtonian expansion of 
h(t) in powers of V is itself an accurate approximation 
to the gravitational-wave signal. 



II. CALCULATION OF THE FOURIER 
TRANSFORM 



We begin with the time-domain signal of Eq. ( |1 . 1| ) , and 
we assume that the signal begins abruptly at a time t m i n 
and ends abruptly at a time t max . We let F min = F(t m i n ) 
and -FJnax = ^(tmax) be the corresponding instantaneous 
frequencies. The relation between F and t is obtained by 
integrating 



dF 
~dt 

which leads to 

t(F) = t c 



!)(> 



57rM [ 



(ttMF) 



11/3 



5M 
256 



(irMF) 



-8/3 



(2.1) 



(2.2) 



where t c ( "time at coalescence" ) is a constant of integra- 
tion. The phase function is then given by 



<!>(F) = $ c -^-(irMF)- 5 /\ 
lb 



(2.3) 



where $ c ( "phase at coalescence" ) is another constant of 
integration. 

The Fourier transform, 



h(f) = / hifye 2 ^ 1 dt, 



(2.4) 

is evaluated by introducing a new integration variable, 
x=(ttMF)^ 3 , (2.5) 



which can be related to t via Eqs. (2.1) and (2.2). After 
some re- arrangement, we obtain 
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where x min = (ttMF^) 1 / 3 , x max = (ttMF^) 1 / 3 , and 



v= (vrM/) 1/3 < 1. 
We have introduced the Fourier integral 



I(v 
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where 



(f>(x;v) 
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It is easy to check that <j)(v;v) = (f>'(v;v) = 0, where a 
prime denotes differentiation with respect to x. The func- 
tion 4>(x;v) initially decreases from 4> m i n (v) = 4>(x m i n ;v) 
to zero as x increases from x m i n to v, and then increases 
from zero to ^ max («) = <^>(a; max ; v) as x increases from v 
to x max . 
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The Fourier integral is evaluated by using (f> as the 
integration variable. Because <f> is not monotonic in the 
interval [a; m i n , a; max ], the integral must be broken down 
into two parts. The first part covers the interval [a; m m, v), 
while the second part covers the interval (v,x max \. It is 
easy to check that the integral can be expressed as 



I{v, x 



\- 16 
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^5 ^min ) + Mv 

3 ^max ) 



where 



and 



Ji(v,x min ) 



(2.10) 
(2.11) 

(2.12) 
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We recall that cf> min = cf>(x lai n]v) and <fi m3X 
are functions of v. 

The coordinate transformation x — > <fi is defined im- 
plicitly by Eq. ( |2.9| ) . It will prove sufficient to invert this 
relation in a neighborhood of x — v, or cj> = 0. The 
following relations are established by Taylor expansion: 



x = v =F 



4V30 7/2 1/2 256 , ; 
15 v 45 
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where the upper sign refers to the interval [a; m ; n , v), while 
the lower sign refers to (v, x max ]. We also have 



3/2,1/2 
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in the second interval. 

To evaluate the integrals Ji(v,x m i n ) and J2(^2anax), 
we let (f> — a — i/3 and deform the contour of integration 
into the complex plane. While the original contour is 
along the a axis, from to m in or </> max , we take the 
new contour to be the union of C and C, where C is 
the curve a = with (3 running from to oo, while C 
is the curve a — </> m in or max with (3 running from oo 
back to 0. The contour is completed by joining C and 
C with the curve (3 — oo, with a running from to <^ m ; n 
or max ; this part of the contour does not contribute to 



the integral. The advantage of this choice of contour is 
that the integrand is exponentially suppressed away from 
f3 = 0, ensuring a rapid convergence of the integral. 

We evaluate th e con tributi on fro m C to J\(v, £ m in) by 
substituting Eq. (2.14) into (2.11), replacing <p by —if3 
and using (3 = and (3 — oo as limits. The integrations 
give rise to r-functions, and we obtain 
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A similar calculation also reveals 
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The contribution from C to Ji(v,x m i n ) is calculated 
by letting <f> — (f> m i n — and expressing the function 
f((j>) = x 2 j (v 3 — x 3 ) as a Taylor series about 4>min- Thus, 

/(</>min - j/?) = /(^min) - i/'(0mi n )/3 H is substituted 

into Eq. (2.11), whose limits are replaced by (3 — oo and 
[3 = 0. The resulting integrations are again elementary, 
and we obtain 



where 

g(v,x) 



I6i x 10 (x 3 + 2v 3 ) 

5 I q t3 rp3 I 



A similar calculation also reveals 
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ie 



-l4>n 



K g(v,x max ). 
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(2.19) 



(2.20) 



(2.14) We note that when x <C 1, the function g(v,x) can be 



well approximated by its first term, 



except 



when x ~ v. In this situation, the expansion for g(v,x) 
does not converge, and our method of calculation breaks 
down. Thus, our expression for the Fourier transform 
will be accurate only when v 3 = wM.f is not too close 
to either x min 3 = nMF min or x max 3 = nMF macK . In 
other words, our expression will be inaccurate near the 
boundaries / = F min and / = F max 

It is easy to see from Eqs. ( 2.18| ) and ( |2.19 ) that the 
contribution from C to Ji{v, x m - ln ) vanishes in the limit 
^min ~~ * 0. Similarly, it can be shown that the contri- 
bution from C to J2(v,:Emax) vanishes (as x max -6 ) in 
the limit x max — > oo. [This behavior is not revealed by 
Eq. ( ^.19] ), which gives g(v,x) as a series expansion for 
small values of x. An alternative expression for g(v,x), 
appropriate for large values of x, can easily be obtained.] 
This allows us to conclude that if the boundaries are 
pushed to F mm = and F max = oo, then C" no longer 
contributes to Ji(v,x mm ) and J2(w,a; max ). 
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Gathering the results, Eqs. ( |2.10D and ( |2.16| )-( ^20| ), 
we find that the Fourier integral can be expressed as 



III. COMPARISON WITH DISCRETE FOURIER 
TRANSFORM 
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where 
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Substituting this i nto Eq. (2.6), treating R and the 0(v 5 ) 
term in Eq. (2.21) as small quantities, we arrive at the 
following expression for the Fourier transform: 



The preceding analysis reveals that apart from win- 
dowing issues, the stationary-phase approximation to the 
Fourier transform is extremely accurate: Apart from the 

. (2-21) 

modulations 5A W and 6?p w , h spa (f) differs from h(f) only 
by a small phase drift Sip of relative order v 10 . In this 
section, we firm up this conclusion by comparing /i spa (/) 
to hgt(f), the discrete Fourier transform of the function 
h(t). 

The discrete Fourier transform is evaluated by Fast 
Fourier Transform (FFT), using the routines of Numeri- 
cal Recipes p^ ]. The time series is prepare d a s follows. 
We begin with h(t) as displayed in Eq. 



^ min g(v,:c m i n ) 
(2.22) 



Hf) = fts Pa (/)(i + 5A w yw>+ s ^l 
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with 
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is the stationary-phase approximation to the Fourier 
transform. Our calculation reveals the existence of two 
types of correction terms. The first is 



Q9 

5tP(v) = -v 5 + 0(v w ), 



(2.26) 



which represents a small, but steadily growing phase 
drift. Notice that Sip is of order 0(v 10 ) relative to ip. 
The other correction terms come as a consequence of the 
abrupt cutoffs imposed at F = F min and F = F max . 
They are 



SA w (v) 
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Notice that SA W represents an amplitude modulation, 
while S(f> w is a phase modulation; both oscillate as a func- 
tion of frequency. The suffix "w" indicates that these 
corrections are associated with "windowing" . 



irrelevant factor QM./r set to unity. Thus, 
h(t) = {irMF) 2 ^e- l,s> , 



1.1[ ), with the 



(3.1) 



where F(t) and &(t) are given by Eqs. ( |E2| ) and (pTq). 
This function is assumed to be nonzero only in the in- 
terval F m j n < F(t) < F max . The duration of the signal 
is 



F t(F maK ) ^(-fmin); 

while the total number of wave cycles is 

M= ^[$(F ma x) -$(F mln )]. 



(3.2) 



(3.3) 



The values of h(t) at the endpoints t(F m i n ) and t(F max ) 
do not agree. This is a potential difficulty for the FFT, 
which considers the signal to be periodic with period T. 
To remedy this, we prepare our time series by padding 
h(t) with zeros on both sides. More precisely, we let h(t) 
be zero in the interval < t < F, be equal to the expres- 
sion ([O]) in the interval F < t < 2T, and be zero again 
in the interval 2T < t < 4T. Thus, the effective duration 
of the time series is four times the duration of the actual 
signal. We choose the value of the parameters t c and <i> c 
such that t(F m i u ) = F and $(F' m in) = 0. This particular 
padding of the time series is somewhat ad hoc (it has not 
been carefully chosen to be the smallest needed to make 
negligible the circular correlations of the FFT), but as we 
shall see the final results are very good and optimization 
would be redundant. 

The zero-padded function h(t) is discretely sampled 
at times tk — kAt, where k = 0,1,..., AN — 1 and 
At = 4T/(47V - 1), with 47V denoting the number of 
sampled points. The Nyquist frequency Q is given 
by iNy = l/(2At), and N must be adjusted so that 
/n y > -F max . The FFT returns the Fourier transform 
of the time series, discretely sampled in the frequency 
domain. The frequency resolution is A/ = (47V At) -1 . 

We denote the discrete Fourier transform of the zero- 
padded time series h(t) by hg t (f), and we wish to com- 
pare this to /i sp a (/), the st ation ary-phase approximation 
given by Eqs. ( [2.24 ) and ( 2.25| ). (It is understood that 
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FIG. 1. The solid curve labeled "SPA" is a plot of 
mod[/i spa (/)], and the dashed curve labeled "FFT" is a plot 
of mod[/ifft(/)]. The inset shows the same curves in a smaller 
frequency interval. 
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FIG. 2. The solid curve labeled "analytical" is a plot of 
1 + <L4 W (/), and the dashed curve labeled "numerical" is a 
plot of A Ic i(f). The inset shows the same curves in a smaller 
frequency interval. 



the correct values for t c and $ c are substituted in these 
equations.) To do this we define the relative amplitude 
A Te \ and relative phase ip Te \ by 



and 



Ad{f) =mod h{£t(f)/h spa (f) 



^rei(f) = arg h S t(f)/h spa (f) 



(3.4) 



(3.5) 



where mod(z) = r is the modulus of the complex number 
while arg(z) = 9 is its argument. If hfn(f) and 



i8 



spa(/) were in perfect agreement, then A rc \ = 1 and 



z 
k 

1pr<A = 0. 

Figure 1 shows plots of mod(/ifft) and mod(/i spa ) as 
functions of / for a signal prepared with M. = 1.25 M© 



in the interval between F„ 



40 Hz and F n 



60 Hz. 



(Such a narrow band is not at all typical for inspiral sig- 
nals sought by any interferometer; we use it merely to 
exaggerate the errors caused by the stationary phase ap- 
proximation to the level where they are visible.) The 
duration of such a signal is T — 15.8 s, and the to- 
tal number of wave cycles is J\f = 749. The FFT was 
taken with N — 2048, giving a Nyquist frequency of 
/Ny = 65 Hz. The figure shows that the agreement is 
not perfect: While the stationary-phase approximation 
seems to give the mean curve, the discrete Fourier trans- 
form displays oscillations about the mean, and these grow 
near the boundary points, / = F m i n and / = F max . We 
shall argue that the discrepancy is entirely caused by 
windowing. 

In Fig. 2 we show plots of A re i(/) — defined in Eq. (3.4) 
- and l + SA w (f) — defined in Eq. (fL"7|), We recall that 
8A w (f) represents the amplitude modulation induced by 
windowing. The near-perfect agreement between A TC \(f) 
and 1 + SA w (f) shows that any discrepancy between 



the discrete Fourier transform and the stationary-phase 
approximation must be attributed to windowing. This 
conclusion is confirmed b y F ig. 3, which shows plots of 
V'reK/) — defined in Eq. ( |3.5|) — and 6ip w (f) — defined 
in Eq. (1.9). Again we see a near-perfect agreement, indi- 
cating that windowing accounts fully for any discrepancy 
between h SD Jf) and hg t (f). 



Equations (1.7) and (1.9) give approximate expressions 
for 5A w (f) and 5ip w (f), and we should expect that in 
some situations, there could be noticeable differences be- 
tween these quantities and the numerically-determined 
A le \(f) and ^ re i(/). Figure 4 indicates that such is indeed 
the case when the frequency interval is expanded. Here, 
the signal is prepared with the same chirp mass as before, 
but the frequency interval is now between F m i n = 40 Hz 
and -F m ax = 1300 Hz; other relevant quantities are listed 
in Table 1 below. Although the agreement is no longer 
near-perfect, it is still remarkably good, and this re- 
enforces our claim that any discrepancy between the dis- 
crete Fourier transform and the stationary-phase approx- 
imation is entirely an artifact of windowing. We have ver- 
ified that the intrinsic correction to the stationary-phase 
approximation, Stp(f) given by Eq. (1.6), is irrelevant in 
this frequency interval: This phase shift is just too small 
to be noticeable at these frequencies. 



IV. OVERLAP INTEGRAL 

We have shown that any discrepancy between /ifft(/) 
and h spa (f) can be fully attributed to windowing, which 
induces amplitude and phase modulations in the Fourier 
transform. In this section we show that these modula- 
tions have no significant effect on operations associated 
with matched filtering. 
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I (Hz) 

FIG. 3. The solid curve labeled "analytical" is a plot of 
5ipv,(f), and the dashed curve labeled "numerical" is a plot 
of ipiei(f)- The inset shows the same curves in a smaller fre- 
quency interval. Notice that the error in the phase does not 
accumulate, and that it is always much smaller than 2n. 




I (Hz) 



FIG. 4. The thin, solid curve labeled "analytical" is a plot 
of 5ipv,(f), and the thick, dashed curve labeled "numerical" is 
a plot of tp re \(f). Notice that here also, the error in the phase 
does not accumulate, and is always much smaller than 2tt. 



The standard theory of matched filtering predicts that 
the loss of signal-to-noise ratio incurred when filtering a 
signal hm{f) with a filter /i spa (/) is equal to Q 



O 



(frfftl^spa) 



/ifftl^fft) (frspa|^spa) 



(4.1) 



We will refer to this quantity as the overlap between the 
two expressions for the signal's Fourier transform. An 
overlap close to unity indicates that the filter is an ac- 
curate representation of the signal, and that using this 
filter in analyzing the data will return the largest possible 
signal-to-noise ratio. We use the notation 



(o|6) 



a*(f)kf)+a(f)b*(f) 

Sn(f) 



df, 



(4.2) 



where S n (f) is the spectral density of the detector noise. 

Notice that the point of view expressed here is that 
hfit(f) is an exact representation of the signal's Fourier 
transform, while ft. spa (/) is an approximate filter. How- 
ever, because of its symmetry in these quantities, O also 
represents the loss in signal-to-noise ratio incurred when 
filtering a signal h spa (f) with a filter hst(f). This is the 
opposite point view, in which the stationary-phase ap- 
proximation is viewed as an exact representation of the 
Fourier transform. 

We evaluate the integrals in Eq. ( [4.1[ ) by turning 
them into discrete sums, using the sampled frequencies 
fk = kAf returned by the FFT. Thus, / a(f) df « 
^ fc a(/fc)A/. It is sufficient for our purposes to use a 
simple analytic model for the noise's spectral density. We 
choose a noise curve that roughly mimics the expected 
noise spectrum of the initial LIGO detector, and set |§ 



S„(f)=S (/ //) 4 + 2 + 2(/// ) : 



(4.3) 



for / > 40 Hz, with f = 200 Hz. The value of S 
is irrelevant for our purposes, and S n (f) is taken to be 
infinite below 40 Hz. 

The overlap O is calculated for a number of chirp 
masses. The results are displayed in Table 1. The conclu- 
sion is clear: For all cases, O > 0.997, indicating that the 
amplitude and phase modulations have very little effect 
on matched filtering. In view of the fact that the modula- 
tions oscillate and never get large, especially away from 
F m i n and F mstx where the instrument is most sensitive, 
this is the expected conclusion. 

It should be noted that for most of the binaries listed 
in Table I, the adopted value for -F max exceeds the fre- 
quency at which the last stable orbit is expected to be 
found. (For equal-mass systems, this frequency is given 
approximately by 1910(Mq/ A4) Hz.) Our gravitational- 
wave signals are therefore not realistic at high frequen- 
cies, something which is already made clear by the fact 
that we do not incorporate post-Newtonian corrections 
into our waveforms. Since our purpose here is simply to 
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TABLE I. The last column gives the overlap O between 
stationary phase and FFT waveforms prepared with (identi- 
cal) chirp mass M given in solar masses in the first column, 
initial frequency i 7 !^ given in Hz in the second column, and 
final frequency _F max given in Hz in the third column. The 
waveforms have a duration T given in seconds in the fourth 
column and a number of wave cycles TV given in the fifth col- 
umn. The number of sampled times is 4JV, where N is given 
in the sixth column, corresponding to a Nyquist frequency 
/Ny given in Hz in the seventh column. In all cases the value 
of -F max is sufficiently large that the contribution to the over- 
lap from higher frequencies can be neglected (assuming the 
initial LIGO noise spectrum given in the text). 



M 


-fin in 


Fm ax 


T 




N 


/Ny 


O 




Hz 


Hz 


s 






Hz 




1.00 


40 


900 


34.6 


2200 


2 16 


948 


0.9999 


1.25 


40 


1300 


23.8 


1521 


2 16 


1375 


0.9998 


1.50 


40 


1300 


17.6 


1122 


2 16 


1863 


0.9997 


1.75 


40 


1200 


13.6 


868 


2 1B 


1204 


0.9997 


2.00 


40 


1500 


10.9 


695 


2 15 


1505 


0.9996 


2.25 


40 


900 


8.9 


570 


2 14 


916 


0.9994 


2.50 


40 


1000 


7.5 


478 


2 14 


1091 


0.9994 


2.75 


40 


1200 


6.4 


409 


2 14 


1279 


0.9994 


3.00 


40 


1400 


5.5 


354 


2 14 


1479 


0.9994 


10.00 


40 


1300 


0.7 


48 


2 11 


1374 


0.9972 



establish the validity of the stationary-phase approxima- 
tion, this lack of realism is not too important. To pro- 
duce a realistic waveform at high frequencies, in a regime 
where the slow inspiral gives way to a rapid merger of 
the two stars, is still an open and challenging problem in 
gravitational-wave research. 
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